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Abstract 

We study inflationary cosmology and the late-time accelerated expansion of the universe in 
non-minimal Yang-Mills (YM) theory, in which the YM field couples to a function of the scalar 
curvature. It is shown that power-law inflation can be realized due to the non-minimal YM field- 
gravitational coupling which maybe caused by quantum corrections. Moreover, it is demonstrated 
that both inflation and the late-time accelerated expansion of the universe can be realized in a 
modified YM-F(R) gravity which is consistent with solar system tests. Furthermore, it is shown 
that this result can be realized also in a non-minimal vector-i ? (i?) gravity. In addition, we consider 
the duality of the non-minimal electromagnetic theory and that of the non-minimal YM theory, 
and also discuss the cosmological reconstruction of the YM theory. 
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I. INTRODUCTION 



Recent observations confirmed that there existed the inflationary stage in the early uni- 
verse, and that also at present the expansion of the universe is accelerating [l], 0]. Various 
scenarios for the late-time acceleration in the expansion of the universe has been proposed. 
In fact, however, the cosmic acceleration mechanism is not well understood yet (for recent 
reviews, see 

There exists two approaches to account for the late-time acceleration of the universe. One 
is dark energy, i.e., general relativistic approach. The other is dark gravity, i.e., modified 
gravity approach. Among the latter approaches studied so far, the modifications to the 
Einstein-Hilbert action, e.g., the addition of an arbitrary function of the scalar curvature to 
it, is one of the most promising latter approaches (for a review, see Q). Such a modified 
theory must pass cosmological bounds and solar system tests because it is considered as an 
alternative gravitational theory. 

A very realistic modified gravitational theory that evade solar-system tests has recently 
been proposed by Hu and Sawicki [s| (for related studies, see [il). In this theory, an effec- 
tive epoch described by the cold dark matter model with cosmological constant (ACDM), 
which accounts for high-precision observational data, is realized as in general relativity with 
cosmological constant (for a review of observational data confronted with modified gravity, 



sec 



10|). This theory can successfully explain the late-time acceleration of the universe. 
In Ref. [ 8|| , however, the possibility of the realization of inflation has not been discussed. 
In Refs. | 111 1 121 . |X3| ] . therefore, modified gravities in which both inflation and the late-time 
acceleration of the universe can be realized, following the previous inflation-acceleration 
unification proposal [3], have been presented and investigated. The classification of viable 
F(R) gravities has also been suggested in Ref. [13]. Here, F(R) is an arbitrary function of 
the scalar curvature R. 

Furthermore, there exists another gravitational source of inflation and the late-time 
acceleration of the universe: a coupling between the scalar curvature and matter La- 



grangian [15J, ll6| (see also (17|). Such a coupling may be applied for the realization of 
the dynamical cancellation of cosmological constant 18]. In Refs. fl9l. I20I. 2l| . the criteria 
for the viability of such theories have been considered. As a simple case, a coupling between 
a function of the scalar curvature and the kinetic term of a massless scalar field in a viable 



modified gravity has been considered [22 



Recently, inflation and the late-time acceleration of the universe in non-minimal electro- 
magnetism, in which the electromagnetic field couples to a function of the scalar curvature, 
have been studied in Ref. 23[ by using the analyzing procedure in the electromagnetic field 



considered in Ref. [24|. It is known that the coupling between the scalar curvature and the 
Lagrangian of the electromagnetic field arises in curved s pac etime due to one-loop vacuum- 
polarization effects in Quantum Electrodynamics (QED) [251 ] . As a result, it has been shown 
that power-law inflation can be realized due to the non-minimal gravitational coupling of 
the electromagnetic field, and that large-scale magnetic fields can be generated due to the 
breaking of the conformal invariance of the electromagnetic field through its non-minimal 



gravitational coupling 1 (see also 27|). The mechanism of inflation in this model is as follows. 



In the very early universe before inflation, electromagnetic quantum fluctuations are gener- 



1 In Ref. [261 ] . gravitational-electromagnetic inflation from a 5-dimensional vacuum state has been consid- 
ered. 
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ated due to the breaking of the conformal invariance of the electromagnetic field and they 
act as a source for inflation. Furthermore, also during inflation electromagnetic quantum 
fluctuations are newly generated and the scale is stretched due to inflation, so that the scale 
can be larger than the Hubble horizon at that time, and they lead to the large-scale magnetic 
fields observed in galaxies and clusters of galaxies. This idea is based on the assumption 
that a given mode is excited quantum mechanically while it is subhorizon sized and then as 



it crosses outside the horizon "freezes in" as a classical fluctuation 28]. These large-scale 
magnetic fields can be the origin of the large-scale magnetic fields with the field strength 
10 _7 -10 _6 G on lOkpc-lMpc scale observed in clusters of galaxies 29J (for reviews of cosmic 



magnetic fields, see [30]). Furthermore, it has been demonstrated that both inflation and the 
late-time acceleration of the universe can be realized in a modified Maxwell- _F(i?) gravity 



proposed in Ref. [13j which is consistent with solar system tests. 

In the present paper, we consider inflationary cosmology and the late-time accelerated 
expansion of the universe in non-minimal non-Abelian gauge theory, called the Yang-Mills 
(YM) theory, in which the non-Abelian gauge field (the YM field) couples to a function of the 
scalar curvature, in order to investigate the cosmological consequences of the non-minimal 
gravitational coupling of the YM filed. Furthermore, we consider a non-minimal vector- F(R) 



gravity. In the past studies, inflation driven by a vector filed has been discussed [3l|, |32 |. 
Moreover, as a candidate for dark energy, the effective YM condensate [33, the Born- 
Infeld quantum condensate [35| and a vector field (36. 37. 38. 39. 40] have been proposed. In 
particular, the possibility that the accelerated expansion of the universe is driven by a field 
with an anisotropic equation of state has been considered in Ref. 140]. As a result, we show 



that power-law inflation can be realized due to the non-minimal gravitational coupling of 
the YM field 2 . Moreover, we demonstrate that both inflation and the late-time accelerated 
expansion of the universe can be realized in a modified Yang-Mills- -F(i?) gravity which is 
consistent with solar system tests. Furthermore, we show that this result can be realized 
also in a non-minimal vector-F(i?) gravity. In addition, we consider the duality of the non- 
minimal electromagnetic theory and that of the non-minimal YM theory, and also discuss 
the reconstruction of the YM theory. 

There are several motivations to study non-minimal YM theory. First of all, we show that 
the appearance of such non-minimal terms in the early universe is compatible with current 
formulations of YM theory due to specific choice of non-minimal function. Second, some 
string compactification may lead to effective scalars-YM-Einstein theory (plus higher order 
corrections). In some cases, one can delete scalars in such a way, that extra curvature terms 
(non-minimal ones) appear in front of YM Lagrangian. Third, the celebrated asymptotic 
freedom phenomenon maybe understood as appearance of non-minimal terms at the early 
universe. 

This paper is organized as follows. In Sec. II we consider a non-minimal gravitational 
coupling of the SU (N) YM field in general relativity. First, we describe our model and derive 
equations of motion from it. Next, we analyze the gravitational field equation, and then show 
that power-law inflation can be realized. In Sec. Ill we consider a non-minimal gravitational 



coupling of the SU(N) YM field in a modified gravitational theory proposed in Ref. [13 
We show that in this theory both inflation and the late-time acceleration of the universe can 
be realized. In Sec. IV we consider a non-minimal vector-F(i?) gravity. Furthermore, in Sec. 



2 In Ref. [4l[, the spontaneous generation of chromomagnetic fields at high temperature has been investi- 
gated. 
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V we consider the duality of the non-minimal electromagnetic theory and that of the non- 
minimal YM theory. In addition, in Sec. VI we discuss the reconstruction of the YM theory. 
Finally, summary is given in Sec. VII. We use units in which fee = c = h = 1 and denote the 
gravitational constant 8ttG by k 2 , so that k 2 = 8vr/M P1 2 , where M P1 = G~ 1/2 = 1.2xl0 19 GeV 
is the Planck mass. Moreover, in terms of electromagnetism we adopt Heaviside-Lorentz 
units. 



II. INFLATION IN GENERAL RELATIVITY 

In this section, following the discussion given in Ref . j23[ . we first consider a non-minimal 
gravitational coupling of the YM field in general relativity. 

A. Model 



We consider the following model action: 



Sgr 
£eh 



J d^xyf^g [ £ E h + £ym 



2k 



■ 2 



R. 



-±I(R)F« u F a ^ 



1 + bg 2 In 



[\/2)F- v F^ 



with 



I(R) = 1 
6=1 1 



f(R) 
n 



4 8tt 2 3 



N. 



FL = d,Al - d v Al + f abc A b Al 



(2.1) 
(2.2) 
(2.3) 

(2.4) 
(2.5) 
(2.6) 



where g is the determinant of the metric tensor g^, R is the scalar curvature arising from the 
spacetime metric tensor g^, and £eh is the Einstein-Hilbert action. Moreover, £ ym with 
I(R) = 1 is the effective Lagrangian of the SU(N) YM theory up to one-loop order 42|, 43 



f{R) is an arbitrary function of R, b is the asymptotic freedom constant, F* is the field 
strength tensor, is the SU(N) YM field with the internal symmetry index a (Roman 
indices, a, b, c, run over 1,2, ... , N 2 — 1, and in F^ u F ailv the summation in terms of the 
index a is also made), and f abc is a set of numbers called structure constants and completely 
antisymmetric [Hj]. Furthermore, fi is the mass scale of the renormalization point, and a 
field-strength-dependent running coupling constant is given by (iij 



1 + bg 2 In \X/n 4 



(2.7) 



where 



(2.8) 



4 



Hence, g is the value of the running coupling constant when X = fi A . 

The field equations can be derived by taking variations of the action in Eq. (12. ip with 
respect to the metric g^ v and the SU (N) YM field A ^ as follows: 



with 



T^ M) = I(R) [g aP F; p F a va e - \g, u T 



(2.9) 



+ - {f'(R)TR, u + g, u D [f'(R)T] - V M V, [f'(R)F]} 



1 + bg 2 In 



JZ — pa pa\xv 



(1/2) F%F^ 



l+bg 2 \n 



X 
7 



1 + bg 2 In 



(1/2) F" 



2X ( 1 + bg 2 In 



X 



and 



=0„ [^gI{R)eF a ^] - I{R)ef abc A b l F c ^ = . 



(2.10) 
(2.11) 

(2.12) 
(2.13) 



where the prime denotes differentiation with respect to R, V M is the covariant derivative 
operator associated with g^, and □ = g^ u V M Vj, is the covariant d'Alembertian for a scalar 

field. In addition, R^ v is the Ricci curvature tensor, while T^ M ^ is the contribution to 
the energy- momentum tensor from the SU(N) YM field. Moreover, e is a field-strength- 
dependent effective dielectric constant (43[, and e ~ 2.72 is the Napierian number. In 
deriving the second equalities in Eqs. (12.111) and (12.121) . we have used Eq. (12.81) . 

We assume the spatially flat Friedmann-Robertson- Walker (FRW) spacetime with the 
metric 



ds 2 



-dt 2 + a 2 (t)dx 2 = a 2 (r])(-dr] 2 + dx 2 ) , 



(2.14) 



where a is the scale factor, and 77 is the conformal time. In this spacetime, g^ v 
diag (— 1, a 2 (t), a 2 (t), a 2 (t)), and the components of R^ and R are given by 



R 



no 



-3 [H + H 2 , Rot = , R^ = [H + 3H Z gij , R = 6 [H + 2H 



(2.15) 



where H = a/a is the Hubble parameter. Here, a dot denotes a time derivative, ' = d/dt. 



B. Power-law inflation 



The (//, v) = (0,0) component and the trace part of the (fi, v) = component of 

Eq. (12.91) . where % and j run from 1 to 3, read 

.2 ( r / \ / \ 1 

T 



H 2 + J 1 = y ^I(R) (bg 2 X + eY) + 3 [-f'(R) (H + H 2 ) + Qf"(R)H (# + AHH 
+ 3f'(R)HA , 



(2.16) 



J, = l -F{R) - F'(R) I 



(2.17) 



and 



2H + 3H 2 + J 2 = I(R)X ( ~e + bg 



~2 



f'(R) [H + 3H' 



+ 6f"(R) [H + 7HH + AH 2 + 12H 2 H ) + 36 f"(R) [H + AHH 



+ 3 



f'(R)H + Af"(R) (h + AHH\ T + f{R)T 



(2.18) 



J 2 = -F'{R)(H + 3H 2 ) + 6F"(i2) 



h + a(h 2 + hh) 



+ 36F"'{R) (h + AHhY 



respectively. Here, X and Y are given by 

X = \E a( - pTopcr \t)\ 2 — |_B a ( propcr )(t)| 2 

Y = |£ ,a (P roper )(t)| 2 -)_ j^a(propcr) |2 



(2.19) 

(2.20) 
(2.21) 



where i^ a (P ropcr )(£) and Bf^ pTopcr \t) are the quantities corresponding to proper electric and 
magnetic fields in the SU(N) YM theory, respectively. In this paper, because we consider 
the case in which there exist the YM electric and magnetic fields as background quantities 
at the 0th order, we here consider that the YM electric and magnetic fields do not have 
the dependence on the space components x. Moreover, Ji and J% are correction terms in a 
modified gravitational theory described by the action in Eq. (13.11) in the next section. Hence, 
because in this section we consider general relativity, i.e., the case F(R) = in the action 
in Eq. (13.21) . here both J\ and J2 are zero. Furthermore, in deriving Eqs. (I2.16P and (12.181) . 
we have used equations in (I2.15P and the following equation: 



F=-2X(e-b?), 



(2.22) 



which follows from Eqs. (J22D, fl2TTTj) and f[2TT2j) . 

In the search of exact solutions for non-minimal YM (electromagnetic)-gravity theory 



sec 



45|, |46|), the problem of off-diagonal components of YM (electromagnetic) stress tensor 
being non-zero while the right-hand side of Einstein equations is zero. In our case, we 
consider as follows. As a simple case, we can consider the following case in which the 
off-diagonal components of T^ M ^ in Eq. (I2.10p vanishes: (i) Only (YM) magnetic fields 
are generated and hence (YM) electric fields are negligible, (ii) B a = (B^B^B^), where 
= B 2 = 0,-Bg 7^ 0, namely, we consider the case in which only one component of 
B a is non-zero and hence other two components are zero. In such a case, it follows from 
divi? a = that the off-diagonal components of the last term on the right-hand side of 



i.e., V M V„ [f'{R)T} are zero. Thus, all of the off-diagonal components of T^ M ^ are zero. In 
this paper (including Sees. Ill and IV) we consider the above case in order to investigate 
the cosmological consequences of the non-minimal gravitational coupling of the YM filed. 
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In Eq. ( I2.13p . because the amplitude of is small, we can neglect the higher order than 
or equal to the quadratic terms in A a ^ and investigate the linearized equation of Eq. (I2.13P 
in terms of A a . The linearized equation of motion in the Coulomb gauge, d^A^(t,x) = 0, 
and the case of A%(t, x) = 0, reads 



A^, x)+[H+j \ A*{t, x) - 1 ( A A?(t, x) 



(2.23) 



(3) 



where A = d l di is the flat 3-dimensional Laplacian. It follows from Eq. (12.231) that the 
Fourier mode A%(k, t) satisfies the equation 







(2.24) 



(2.25) 



I I a' 

Replacing the independent variable t by rj, we find that Eq. (I2.24p becomes 
d 2 Af(k, V ) 1 dI{ri) dAt{k,rj) 2 

By using the WKB approximation on subhorizon scales and the long-wavelength approxi- 
mation on superhorizon scales, and matching these solutions at the horizon crossing 24| , we 
find 

2 



IA a (M)| 2 = |C(fc)p 



i 



2kI( Vk 



1 



1 1 



dl(rjk 



% lyn) 



(2.26) 



2kI( Vk ) drj 

where rjk and r]t are the conformal time at the horizon-crossing and one at the end of 
inflation, respectively. Consequently, from Eq. ( 12.261) we obtain the amplitude of the proper 
YM magnetic fields on a comoving scale L = 2ir/k in the position space 

k\C(k)\ 2 k 4 



| jga (proper) I 



7P 



1 + I r bc u b u ck\C{^l 
2 J 2n 2 



(2.27) 



where u b (= 1) and u c (= 1) are the quantities denoting the dependence on the indices b 
and c, respectively. Thus, from Eq. (12.271) we see that the YM magnetic fields evolves as 
[^(proper) ^ 1 2 = \B a \ 2 /a\ where \B a \ is a constant. 

In this case, we find that Eqs. (12.161) and (I2.18P are reduced to 



H 2 = K 2 \^-I{R){e-b~g 2 " 



f'(R) (H + H 2 ) + 6f"(R)H (ii + AHh)] (e - bg 2 ) 



- Af'(R)H 2 e 



\B 



a\2 



(2.28) 



and 



1 



f\R) [H + 3H' 



2H + 3H 2 = k 2 1 -I(R)(e - 3bg 2 ) 
6 



+ Gf"(R) [H + 7HH + AH 2 + 12H 2 H) + 36/' "(R) [H + AHH 



[e - bg 



+ 4 



f'(R) 



■H + H' 



12 f "' (R)H (H + AH H J e + lQf'(R)H 2 bg 



\B 



a\2 



(2.29) 



respectively. Eliminating I(R) from Eqs. (12.281) and (I2.29p . we obtain 



H + t^w r2 = k2 I f{R) i ~ {2e + h ~ g2)H 



,~,2 



3e: — 75g 



7.2 



e + 8&# 



e-bg" 

f 3/"(i2) (e - bg 2 )H - 2(e + 2bg 2 )HH + 4(e - &# 2 )# 2 - 24eH 2 H 

> 15012 



+ 18 f"(R)(e - bg 2 ) H + 4M 



(2.30) 



We here note the following point. From Eq. (12.111) . we see that the value of e depends on 
the field strength, in other words, it varies in time. In fact, however, the change in time of e 
is smaller than that of other quantities because the dependence of e on the field strength is 
logarithmic, so that we can approximately regard e as constant in Eq. (I2.30p . (Thus, from 
this point we regard e as constant.) 

Here we consider the case in which f(R) is given by the following form: 



f(R) = f m (R) 



c\ (R/m 2 



c 2 (R/m 2 ) n + 1 



which satisfies the conditions: 



lim fns(R) = — — const 
hmf HS (R) = 0. 

H — n) 



(2.31) 

(2.32) 
(2.33) 



Here, C\ and c 2 are dimensionless constants, n is a positive constant, and m denotes a mass 
scale. This form, fns(R), has been proposed by Hu and Sawicki The second condition 
(I2.33P means that there could exist a flat spacetime solution. Hence, because in the late 
time universe the value of the scalar curvature becomes zero, the YM coupling I becomes 
unity, so that the standard YM theory can be naturally recovered. 

In order to show that power-law inflation can be realized, we consider the case in which 
the scale factor is given by a(t) = a {t/t) p , where i is some fiducial time during inflation, a 
is the value of a(t) at t = i , and p is a positive constant. In this case, H = p/t, H = —p/t 2 , 
H = 2p/t 3 , and H = —Qp/t 4 . Moreover, it follows from the fourth equation in (I2.15P that 
R = 6p(2p — l)/t 2 . At the inflationary stage, because R/m 2 ^> 1, we are able to use the 
following approximate relations: 



hs(R) 



Cl 
(■■■2 



1 

C-2 



R 



77r 



(2.34) 



Substituting the above relations in terms of a, H and R, and the approximate expressions 
of fu S (R), fns( R ) and fus( R ) derived from Eq. (I2T3D into Eq. fTO0|) . we find 



P 



tP 



n + 1 



-Ci) (n + 



l)e 2 + 3(n - l)bg 2 e + 6 



,7,2\2 



3 n+1 n n (n + l) n+1 c 2 2 



(n - 1)5 + 2bg 



a2 



\B a \ 2 K 2 m 2n 



1/1 



(2.35) 



(2.36) 



S 



Hence, if n 3> 1, p becomes much larger than unity, so that power-law inflation can be 
realized. Consequently, it follows form this result that the YM field with a non-minimal 
gravitational coupling in Eq. (12.31) can be a source of inflation. This result is the same as in 
non- minimal Maxwell theory |23| . 

In this paper we consider only the case in which the values of the terms proportional to 
f(R), f"(R) and f'"(R) in the right-hand side of Eqs. (JOED and are dominant to 

the value of the term proportional to I(R). Among the terms proportional to f'{R), f"{R) 
and f"'(R), the term proportional to f'{R) is dominant, and its value is order f'(R)H 2 « 
n (C1/C2 2 ) (H 2 /m 2 ) (R/m 2 ) < - n+1 ' ) ; which can be derived by using Eq. (I2.34p . Here, it follows 
from H = p/t and R = Qp(2p — l)/t 2 that R is order 10H 2 . The condition that the 
term proportional to f'{R) is dominant in the source term would be I(R)/ [f'(R)H 2 ] ~ 
10c2(R/m 2 ) n jn <C 1. This would require extremely small C2 because at the inflationary 
stage R/m 2 ^> 1 and n 1. In such a case, the value of the right-hand side of Eq. (12.301) . 
which is order k 2 f (R)H 2 \B a \ 2 j a 4 , can be order H 2 . Consequently, the right-hand side of 
Eq. (I2.30p can balance with the left-hand side of Eq. (I2.30p . and hence Eq. (12.301) can be 
satisfied without contradiction to the result, i.e., power-law inflation in which p is much 
larger than unity can be realized. The reason why we consider the case in which the term 
proportional to I{R) on the right-hand side of Eqs. (I2.28P and (I2.29P is so small in comparison 
with the term proportional to f'{R) that it can be neglected is as follows (23[: If the 
opposite case, namely, the term proportional to I{R) is dominant to the term proportional 
to f'(R), Eqs. (T2T28]) and (I2T29D are approximately written as H 2 « (l/6)K, 2 I(R)\B a \ 2 /a 4 and 
2H + 3H 2 « (l/6)K 2 I(R)\B a \ 2 /a 4 , respectively. Thus, in this case it follows from Eqs. (|2T28|) 
and (12.291) that H 2 and 2H + 3H 2 are the same order and their difference, 2H + 2H 2 , must 
be much smaller than H 2 . In fact, Eq. (I2.30p implies that H + [e/ (e — bg 2 )} H 2 balances 

with much smaller quantity than k 2 I (R)\B a \ 2 / a A . Now, {if + [e/ (e - bg 2 )} if 2 } / H 2 = 

ej (e — bg 2 ) — 1/p <^ 1 and hence p must be smaller than unity because e > and b > 0. 
Consequently, in this case power-law inflation cannot be realized. 

Finally, we note the following point. The constraint on a non-minimal gravitational 
coupling of matter from the observational data of the central temperature of the Sun has been 
proposed [2l[ . Furthermore, the existence of the non-minimal gravitational coupling of the 
electromagnetic field changes the value of the fine structure constant, i.e., the strength of the 
electromagnetic coupling. Hence, the deviation of the non-minimal electromagnetism from 
the ordinary Maxwell theory can be constrained from the observations of radio and optical 
quasar absorption line s |47l| . those of the anisotropy of the cosmic microwave background 
(CMB) radiation [48|, l49j. those of the absorption of CMB radiation at 21 cm hyper fine 
transition of the neutral atomic hydrogen [50| , and big bang nucleosynthesis (BBN) [5lL 52l 
as well as solar-system experiments [53( (for a recent review, see [54j). On the other hand, 
because the energy scale of the YM theory is higher than the electroweak scale, the existence 
the non-minimal gravitational coupling of YM field might influence on models of the grand 
unified theories (GUT). 



III. INFLATION AND LATE-TIME COSMIC ACCELERATION IN MODIFIED 
GRAVITY 



Next, in this section we consider a non-minimal gravitational coupling of the YM field in 
a modified gravitational theory proposed in Ref. 13 . 
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We consider the following model action: 



S 



MG 



d 4 xy/^g [ £mg + £ym 



c 



MG 



1 

2k 



2 [R + F(R)] 



(3-1) 



(3.2) 



where F(R) is an arbitrary function of R. Here, £ym is given by Eq. (12.31) . We note 
that F(R) is the modified part of gravity, and hence F(R) is completely different from the 
non- minimal gravitational coupling of the YM field f(R) in Eq. ( 12.4ft . 

Taking variations of the action Eq. (13.11) with respect to the metric g^, we find that the 
field equation of modified gravity is given by [13| 



[1 + F'(R)} R, u - l -g, v [R + F(R)} + g, v UF\R) - V^ V F'{R) = k 2 T™ . (3.3) 

The (/i, v) = (0,0) component and the trace part of the (//, i>) = component of 

Eq. (13.31) . where % and j run from 1 to 3, are given by Eqs. (12.161) and (12.18p . respectively. 

Here we consider the same case as in the preceding section. In this case, eliminating I(R) 
from Eqs. (12.161) and 02.18p . we obtain 



H + 



s-bg 



^H 2 + 



6(5 — bg 



7,2 



■F(R) - F'(R) 



bg 



7,2 



b~g 



7,2 



H 



b~g 2 



H' 



3F"(R) 



if + 4 H 



HH )\ 



18F"'(R) (H + 4HH\ 



K 2 (f'(R)t-(2e + b~g 2 )H + 



7,2 



3s — 7bg 



~2 



e + Sbg 



e — bg 2 

f 3f"{R) (e - bg 2 )H - 2(e + 2bg 2 )HH + 4(e - b~g 2 )H 2 - 2AeH 2 H 

\ B a\2 



+ 18 f'"(R)(e - bg 2 ) \H + AHH 

Here we consider the case in which F(R) is given by 

[(R/M 2 ) - (R /M 2 )f +1 + (R /M 2 f l+1 



F(R) = -M l 



c 3 + c 4 [[(R/M 2 ) - (R /M 2 )] 2l+1 + (R /M 2 r i+i j 



\2l+l\ 



(3.4) 



(3.5) 



which satisfies the following conditions: lim^oo F(R) = —M 2 jc A = const, liniR^o F(R) = 
0. Here, C3 and C4 are dimensionless constants, / is a positive integer, and M de- 
notes a mass scale. We consider that in the limit R — > 00, i.e., at the very early 
stage of the universe, F(R) becomes an effective cosmological constant, lim^oo F(R) = 
-M 2 /c A = -2Ai, where A s (> H 2 ) is an effective cosmological constant in the very 
early universe, and that at the present time F(R) becomes a small constant, F(Rq) = 



—M 2 {Ro/M 2 f +1 j c 3 + c 4 (R /M 2 ) 2l+1 = -2Rq, where R (« H 2 ) is current curva- 
ture. Here, H is the Hubble constant at the present time: H = 100/i km s _1 Mpc -1 = 
2.1h x 10" 42 GeV « 1.5 x 10" 33 eV 0, where we have used h = 0.70 [56 . 



10 



Furthermore, we consider the case in which f(R) is given by the following form: 

tK ) /N ° l ] " c 5 + c 6 {[(i?/M 2 ) - {R,/M 2 )f^ + (R /M 2 ) 2 ^} ' ^ 

which satisfies the following conditions: lim^oo /no(-R) = V c 6 = const, liniR^o /no(-R) = 0. 
Here, C5 and C6 are dimensionless constants, and q is a positive integer. The form of F(R) 
in Eq. (13.51) and f^o(R) in Eq. (13.61) is taken from Ref. [13] . This form corresponds to the 
extension of the form of fns(R) m Eq. (12.311) . It has been shown in Ref. [l3| that modified 
gravitational theories described by the action (13.21) with F(R) in Eq. (13. 5p successfully pass 
the solar-system tests as well as cosmological bounds and they are free of instabilities. 
Making the same considerations as in Ref. j23|, we find that at the very early stage 

of the universe, it follows from Eq. (13.41) that a(t) oc exp ^A/A;/3tj , so that exponential 

inflation can be realized, and that at the present time, it follows from Eq. (13 .4ft that a(t) oc 

exp (^J R^/St^ , so that the late-time acceleration of the universe can be realized. These 

results are also the same as in non-minimal Maxwell-F(i?) gravity 23[ . 

Finally, we note the following point about the logarithmic contribution to modified grav- 
ity, namely, the case in which the Lagrangian of modified gravity in Eq. ( 13.21) are given by 
£mg = 1/ (2^ 2 ) [R + F(R) + In (R/M 2 )]. Following to the considerations in the previous 
subsections, because the logarithmic term is sub-leading contribution, in also this case both 
inflation and the late-time acceleration of the universe can be realized. The qualitative 
difference from the case of the previous subsections is only that in the limit R — > 00 the 
gravitational modification term, F(R) + ln (R/M 2 ) does not become constant. In fact, how- 
ever, if it is considered that some cut off scale of R in the very early universe should exist, 
the logarithmic contribution does not diverge in this limit, and hence the cosmology of this 
case is the same as that of the previous sections. 

IV. NON-MINIMAL VECTOR MODEL 

In this section, we consider the cosmology in the non-abelian non-minimal vector-F(i?) 
gravity. 

We consider the following model action: 



5mg = J d %V^g [ Cmg + £>v ] , (4.1) 

£v = I(R) {~\f; v F^ - V[A a2 }} , (4.2) 

where C M g is given by Eq. fl3T2|) . F£ v is given by Eq. fl2T6|) . and A a2 = g^A^A*. (As the 
generalization of the above non-minimal vector model, one can consider a model in which 
the derivative in is the gauge covariant derivative given by = — igA^, where 
A^ = A^t 11 . Here, r a are matrices and their commutation relations is conventionally written 
as the standard form [Y a , r b ] = if abc r c j44j. In the present paper, however, as a simple 
non- minimal vector model we consider the theory described by the action in Eq. (14. 2p ). 

We should note that the last term y[v4 a2 ] in the action (14.21) is not gauge invariant but 
can be rewritten in a gauge invariant way. For example if the gauge group is a unitary 
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group, we may introduce a a-model like field U, which satisfies U'U = 1. Then the last 
term could be rewritten in the gauge invariant form: 

V[A a2 } -> V [ctr (rfA°U) (U^A^U)] . (4.3) 

Here c is a constant for the normalization. If we choose the unitary gauge U — 1, the term 
in (14.31) reduces to the original one: y[A a2 ]. This may tells that the action (14.21) described 
the theory where the gauge group is spontaneously broken. 

The field equations can be derived by taking variations of the action in Eq. (14. ip with 
respect to the metric g^ v and the vector field A^ as follows: 

[1 + F'(R)} R„ v - [R + F(R)} + g, u UF\R) - V,V V F'(R) = k 2 T$) } (4 . 4) 



with 



= I{R) [g aP F^ a + 2 ^^1 - \g^ 

+ I {f{R)FR» v + g, v D [f'(R)f] - V M V„ [f'(R)f] } , (4.5) 



T = F a F a » v + AV[A a2 ] , (4.6) 



and 



p^fy [V^gI(R)F a ^ u ] - I(R) ^bc A b^ F c„u + 2 dV d ^p A au ^j = 



(4.7) 



where T$ is the contribution to the energy-momentum tensor from A a . 



Here, as an example, we consider the case in which V[A a2 ] is given by a class of the 
following power-law potential: 

( A a2 \ n 

V[A a2 } = V^—j , (4.8) 

where V is a constant, fa denotes a mass scale, and n{> 1) is a positive integer. 

Similarly to Sec. II B, because the amplitude of A a ^ is small, we neglect the higher order 
than or equal to the quadratic terms in A a ^ and consider the linearized equation of Eq. (14. 7p 
in terms of A a . For the power-law potential given by Eq. (14.81) . the linearized equation of 
motion under the ansatz d 3 Aj(t, x) — and A%(t, x) = is the same as Eq. (I2.23p . 3 

The = (0,0) component and the trace part of the (//, v) = component of 

Eq. (14.41) . where % and j run from 1 to 3, read 



H 2 + qF(R) - F'(R) (k + H 2 

2 / 

= y ( I(R) {Y + 2V[A a2 } } + 3 [-f'(R) (H + # 2 ) + Gf"(R)H (h + 4m) 



1f{R)RT\ , (4.9) 



This is similar to the Coulomb gauge but since the action (|4. 2|) is not gauge invariant, or gauge symmetry 
is completely fixed by the unitary gauge as in after (|4.3[) . this condition is only a working hypothesis. 
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and 



2H + 3H 2 + -F(R) - F'(R) [H + 3H 



6F"{R) 



H + A(H 2 + HH 



3QF"'(R) [H + AHH 



+ 



+ Qf"(R) (li + 7HH + AH 2 + 12# 2 if) + 36f"'(R) [h + AHH 



f'(R) [H + 3H'' 



+ 3 



f'(R)H + Af"(R) (h + 4M)j £ + f'(R)T\ , 



(4.10) 



respectively. In deriving Eqs. ( 14. 9 p and ( 14.101) . we have used equations in ( 12. 15(1 . 

Here we consider the same case as in the previous sections. Moreover, we here consider 
the case in which = 0. In this case, we have (1/a 2 ) AfAfdV[A a2 ]/ (dA a2 ) = nV[A a2 ]. 
Consequently, using this relation, we find that Eqs. (14. 9 p and (14.101) are reduced to 



H 2 + -F(R) - F'(R) [H + H 2 
6 



-I(R) - f'(R) (h + 5# 2 ) + 6f"(R)H (h + AHH 



IB 



a\2 



+ l-I(R) -2f'(R) H + (l + 2n)H 2 + 12 f"(R)H [H + AHH) } V[A 



\a2- 



(4.11) 



and 



2H + 3H 2 + ^F{R) - F\R) (h + 3H 2 } 
+ 6F"(R) 



H 



H + A(H 2 + HH 
1 



6 



f 36 F"'{R) (Ji + AHH^ 
I(R) + f'(R) ( -5H + H 2 ) + 6f"(R) (h-HH + AH 2 - 20H 2 H 





\B a \ 


(h + ahhY 


a 4 



I(R) - 2f\R) (1 + 2n) H + (3 + Qn - An 2 ) H 2 



+ \2f"{R) H + {7- An) HH + AH 2 + A(3- An) H 2 H 



72 f"'(R) [H + AHH) >V[A 



\a2- 



(4.12) 



13 



respectively. Eliminating I(R) from Eqs. (14.111) and (I4.12p . we obtain 



H + H 2 +{ -F(R) - F'(R)H 2 + 3F"{R) 
6 



H + a(h 2 + HH) + 18F"'{R) (H + AHH\ ' 



f{R) ( -2H + 3H 2 ) + 3/"(i2) ( H - 2HH + AH 2 - 2AH 2 H 



\B 



a\2 



2\ -f(R) nH+(l + 2n- 2n 2 ) H 1 



+ 18 f\R) (H + AHH 



+ 3f"(R) H + 2 (3 - 2n) HH + AH 2 + 8(1- 2n) H 2 H 



18 f" '(R) ( H + AHH^j 2 



a2" 



(4.13) 



In the case that |5f (proper) (t)| 2 = |5 a | 2 /a 4 , \/[A a2 ] oc a~ 2n . Hence, if n = 2, the time 
evolution of V[A a2 } is the same as that of \Bf propei \t)\ 2 . On the other hand, if n > 2, 
y^" 2 ] decreases much more rapidly than | J B j a (P r °P cr )(t)| 2 during inflation. Thus, in the latter 
case we can neglect the terms proportional to V[y4 a2 ] on the right-hand side of Eq. (14.131) . 
Consequently, it follows from Eq. (14.131) that when we consider the case in which similarly 
to the preceding section, F(R) and f(R) are given by Eqs. (13. 5p and (13.61) . respectively, we 
can make the same consideration as the preceding section, and hence power-law inflation 
and the late-time acceleration of the universe can be realized. 

Furthermore, as another case, we consider the case in which F(R) is given by Eq. (13.51) 
and f{R) is given by the following form: 



f(R) = RR) 



c 7 (R/M 2 y-1 
c 8 [R/M 2 Y + 1 



(4.14) 



which satisfies the following conditions: lim^oo/(_R) = Cj/cs = const, liniR^o /(-R) = — 1. 
Here, C7 and c§ are dimensionless constants, q is a positive constant, and M denotes a mass 
scale. In this case, the value of I(R) = 1 + f(R) becomes close to zero when that of R is 
very small, namely, at the present time. Making the same consideration as the preceding 
section, we can also find in this case that power-law inflation and the late-time acceleration 
of the universe can be realized. 



V. DUALITY 

In this section, we consider the duality of the non-minimal electromagnetic theory and 
that of the non-minimal YM theory. 



A. Duality of the non-minimal electromagnetic theory 



We consider the duality of the action of /(-R)-coupled electromagnetic theory: 
1 



5 



fA 



d 4 x^f(R)F AlMU F 



/<!/ 

A ' 



A fiu 



= d^Ay - d v A 



(5.1: 
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where f(R) is an arbitrary function of R and A u is the U(l) gauge field. Before going to 
the duality of the action in Eq. (15.11) . we consider the duality without gravity: 

S A = \J d*xF Aflu F% , F Apu = d p A u - d u A^ . (5.2) 

By introducing a new field B^, the action can be rewritten as 



SpB = \j d 4 x + X -e^° [d p B v ) F p )j . 



(5.3) 



Here, F pu is an independent field (not given in terms of A p or B p as in F Apu ). The variation 
of Bfj, gives 

e^ pa d u F pa = , (5.4) 
which tells that F pu can be given in terms of a vector field A p as 

F l „, = F AlJV = d IM A v -d v A li . (5.5) 

Then the action Sp B in Eq. (15.31) reduces to S A in Eq. (15.21) . On the other hand, by the 
variation of F pu , we obtain 



F^ = --e^ pa d p B a . (5.6) 



1 

4 

By substituting Eq. (15.61) into the action in Eq. (15. 2p . we obtain 

Sb = \J d 4 xF Bpu F% , F Bltv ee d p B u - d v B, . (5.7) 
Eqs. (15. 5p and (15. 6p give 

FT = \e^F Bpa , (5.8) 

which tells that Fb^ v are dual to F Apu , that is, the magnetic field exchanges with the electric 
field. 

We now consider the action in Eq. (15.11) . which can be rewritten as 

S fPB = ~ fd A x {j=gf{R)F ia ,F'» + l -e^ (d p B u ) F pCT | . (5.9) 

Now F pu is an independent field again. From the variation of B p , we obtain Eq. ( 15. 4p . which 
can be solved as Eq. (I5.5p . and we find the action in Eq. (15.91) is equivalent to Eq. (15.11) . 
On the other hand, by the variation of F^, instead of Eq. (15.61) . we obtain 

F^ = —„ ^£l. (5.10) 

*f(R)V=g 

Then substituting Eq. (I5.10p into Eq. (15. 9p and using the identity 

e a ^e^ 5 g pa g^ = 2g {g^g° s - g p5 g^) , (5.11) 
we obtain an action dual to Eq. (15. ip : 

SfB = \J ^ x ^J7gj F Bv> F B , F b,u = d p B u - d v B p . (5.12) 
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B. Duality of the non-minimal Yang-Mills theory 

As in case of the electromagnetic theory, we may consider the duality of the action of 
/(-R)-coupled Yang-Mills theory: 

S fA = \j d 4 x^f(R)F Allu F^ , F Aflu = d,A u - d u A, + f abc A b ^A c u . (5.13) 

The action can be rewritten in the following form 

S fFB = \ J d A x {v=-gf{R)F, v F^ + \^F%F^ . (5.14) 

Now F^ is an independent field again. From the variation of A ^ in F® , we obtain 

e^ pa D v F prT = . (5.15) 

Here is a covariant derivative. The solution of (15.151) is given by 

= F AfMV . (5.16) 

By substituting (15.161) into (I5.14p . we obtain (15.131) . On the other hand, by the variation of 
F^u, we obtain 

= —~ . (5.17) 

Then substituting Eq. (15.171) into Eq. (15.141) and using the identity (15. lip , we obtain a dual 
action: 

S f B = \j *x^-±-F Aia ,F7 . (5.18) 
Note that dual form of the action maybe useful in the cosmological considerations. 

VI. RECONSTRUCTION OF THE YM THEORY 

In this section, we indicate how to reconstruct the YM theory from the known universe 
evolution (for a review, see (57|). 

We now consider the following action: 

S = J d*xy/^ (A + t (i^F^)) . (6.1) 

By introducing an auxiliary scalar field 0, we may rewrite the action (16. ip in the following 
form: 

s = J d'xV^ (A + ±p{<p)f; v f^ + \q{4>)} . (6.2) 

By the variation of 0, we obtain 

= P'(<j ) )F^F a ^ + Q\ ( f ) ) , (6-3) 
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which could be solved with respect <p as (p = (ft [F®F atJ,v ^. Here, the prime denotes differ- 
entiation with respect to <p. Then by substituting the expression into (16.21) we obtain the 
action (16. ip with 

T (F^F a ^) = \{P{<f> (F^F^)) F; v F a ^ + Q (0 {F; u F a ^)) } . (6.4) 

By the variation of the action (16. 2p with respect to the metric tensor g^, we obtain the 
Einstein equation: 

± (r, v - l -Rg^ = - l -P(<p)F« p Fy + l -g, v {P(<p)F^F a ^ + Q{<j>)) . (6.5) 

On the other hand, by the variation with respect to A^, we obtain 

= d v (y/=g~P(4>)F a ^ 1 ) - y^g~P(4>)f abc A b u F c ^ . (6.6) 

For simplicity, we only consider the case that the gauge algebra is SU(2), where f abc = e abc , 
and we assume the gauge fields are given in the following form 



Here a is a constant with mass dimension and A is a proper function of t. In general, if 
the vector field is condensed, the rotational invariance of the universe could be broken. In 
case of (16.71) . the direction of the vector field is gauge variant. Then all the gauge invariant 
quantities given by (16. 7p do not break the rotational invariance. 
By the assumption, (I6.3P has the following form: 

= 6 (-a 2 A 2 e 2A a" 2 + a 4 e 4A a" 4 ) P'(0) + Q'(0) , (6.8) 

and (t, t)-component of Eq. (16.51) is given by 

= Atf 2 _ I (a 2 A 2 e 2 V 2 + a W 4 ) - ±Q(<j>) . (6.9) 

The fi = component of (16.61) becomes identity and /i = i component gives 

= d t (aP(0)Ae A ) - 2a 2 a- 1 P(0)e 3A . (6.10) 

Since we can always the scalar field <p properly, we may identify the scalar field with the 
time coordinate <fi = t. Then by differentiating Eq. (I6.9P with respect to t and eliminating 
Q = Q'(4>), we obtain 

= ^HH + a 2 \ 2 e 2X a~ 2 P 

K 

-P [a 2 (AA + A 3 - A 2 #) e 2A a~ 2 + 2a 4 (A - if) e 4A a~ 4 } . (6.11) 

Furthermore by eliminating P by using (I6.10p . we find 

2HH 

P=— ( T . TTT T . (6.12) 

k 2 \2a 2 a~ 2 e 2X (\ 2 + AAJ - a 4 e 4A a" 4 # j 
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(^ = i = 1,2,3) 



Then by using (I6.12p . we can eliminate P (and P) in (16.101) and obtain 



= 2 (AA + A 2 + 3A 2 a) - a'V x a- 2 H + 4 |A 3 + AA — aW 2 //} (A - H) 

+ {2(A' + AA)-«VV 2 /f}{| + | + ff + i + A- 25X^1 . (6.13) 

If we give a proper a = a(t) and therefore H = H(t), Eq. ( 16. 13f) can be regarded as a third 
order differential equation with respect to A. If we find the solution of A with three constants 
of the integration, we find the explicit form of P(<p) = P(t) by using (I6.12p and further obtain 
Q(4>) by using (16.91) . Then we find the explicit form of three parameter families of the action 
( 16. 2p . This tells that almost arbitrary time development of the university could be realized 
by the action Q or (IBTT]) . 

As an example, we may consider the case of the power law expansion: 




(6.14) 



Here t\ and hi are constants. By assuming 

A = (/i 1 -l)ln^+A 1 , (6.15) 

(Ai is a constant), Eq. (16. 3ft reduces to the algebraic equation: 
2h 

= - — X —X 2 + {-Ah\ + 13/ix + 2) X + {hi - if (hi - 2) (AHi - 20) . (6.16) 
hi — 1 

Here X = a 2 t\e . If (I6.16P has a real positive solution with respect to X, we obtain Ai and 
therefore the exact form of A. Then we can reconstruct a model to give the power expansion 
(I6.14p . Similarly, any other universe evolution history maybe reproduced by specific form of 
the action under consideration. 



VII. CONCLUSION 

In the present paper, we have considered inflationary cosmology and the late-time ac- 
celerated expansion of the universe in the YM theory, in which the YM field couples to a 
function of the scalar curvature, in order to investigate the cosmological consequences of 
the non- minimal gravitational coupling of the YM filed. As a result, we have shown that 
power-law inflation can be realized due to the non-minimal gravitational coupling of the YM 
field. Moreover, we have demonstrated that both inflation and the late-time accelerated ex- 



pansion of the universe can be realized in a modified YM-F(R) gravity proposed in Ref. [13 
which is consistent with solar system tests. Furthermore, we have shown that this result can 
be realized also in a non-minimal vector-F(i?) gravity. In addition, we have considered the 
duality of the non-minimal electromagnetic theory and that of the non-minimal YM theory. 
Furthermore, we also discussed the reconstruction of the YM theory from the known uni- 
verse history expansion. As an example, it has been shown that a model to give the power 
expansion of the scale factor can be reconstructed. 
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Finally, we remark the following point. It is interesting that our models maybe extended 
by another gauge-non-invariant non-minimal coupling with the curvature like the ones done 
in Ref. [58[ (for models of vector curvaton, see [59]). Such non-minimal vector curvaton 
may give extra contribution to curvature perturbations if compare with the present models. 
Another important point is related with the exit from the inflation. In the models under 
consideration it maybe realized via the gravitational scenario, as the instability of de Sitter 
universe, due to extra gravitational terms. This scenario will be investigated in detail else- 
where. It maybe also relevant for the study of future universe: if our universe will stay as 
A-CDM one forever or it will evolve to other singular/non-singular state. 

Acknowledgments 

We are grateful to M. Sasaki for very helpful discussion of related problems. The work 
of K.B. was supported in part by the open research center project at Kinki University and 
that by S.D.O. was supported in part by MEC (Spain) projects FIS2006-02842 and PIE2007- 
50/023. This work by S.N. is supported in part by the Ministry of Education, Science, Sports 
and Culture of Japan under grant no. 18549001 and 21st Century COE Program of Nagoya 
University provided by the Japan Society for the Promotion of Science (15COEG01). 



[1] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003); H. V. Peiris 
et al. [WMAP Collaboration], ibid. 148, 213 (2003); D. N. Spergel et al. [WMAP Collabora- 
tion], ibid. 170, 377 (2007). 

[2] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 
(1999); A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998); 
P. Astier et al. [The SNLS Collaboration], Astron. Astrophys. 447, 31 (2006); A. G. Riess et 
al., |arXiv:astro-ph/0611572 

[3] P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003); V. Salmi, AIP Conf. Proc. 
782,166 (2005). 

[4] T. Padmanabhan, Phys. Rept. 380, 235 (2003). 

[5] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006). 
[6] R. Durrer and R. Maartens, arXiv:0711.0077 [astro-ph]. 

[7] S. Nojiri and S. D. Odintsov, Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007) 

[arXiv:hep-th/0601213] ; larXiv:0801 .48431 [astro-ph] . 
[8] W. Hu and I. Sawicki, Phys. Rev. D 76, 064004 (2007). 

[9] S. A. Appleby and R. A. Battye, Phys. Lett. B 654, 7 (2007); S. Nojiri and S. D. Odintsov, 
ibid. 652, 343 (2007) |arXiv:0706. 13781 [hep-th]]; L. Pogosian and A. Silvestri. larXiv:0709.0296l 
[astro-ph]; S. Capozziello and S. Tsujikawa, larXiv: 0712.2268 [astro-ph]; S. A. Appleby and 
R. A. Battye, arXiv:0803.1081 [astro-ph]; S. Tsujikawa, K. Uddin and R. Tavakol, Phys. Rev. 
D 77, 043007 (2008); S. Tsujikawa, K. Uddin, S. Mizuno, R. Tavakol and J. Yokoyama, 
larXiv: 0803 .11061 [astro-ph]. 

[10] S. Capozziello and M. Francaviglia. larXiv:0706.1146l [astro-ph]. 

[11] S. Nojiri and S. D. Odintsov. larXiv:0710.1738l [hep-th], 

[12] G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov, L. Sebastiani and S. Zerbini, 
larXiv:0712.40T7l [hep-th]. 



19 



[13] S. Nojiri and S. D. Odintsov, Phys. Lett. B 657, 238 (2007) |arXiv:0707.194T1 [hep-th]]. 

[14] S. Nojiri and S. D. Odintsov, Phys. Rev. D 68, 123512 (2003) [arXiv:hep-th/0307288| . 

[15] S. Nojiri and S. D. Odintsov, Phys. Lett. B 599, 137 (2004) |arXiv:astro-ph/0403622| ; PoS 

WC2004, 024 (2004) |arXiv:hep-th/0412030| . 
[16] G. Allemandi, A. Borowiec, M. Francaviglia and S. D. Odintsov, Phys. Rev. D 72, 063505 

(2005) [arXiv:gr-qc/0504057] . 
[17] N. Deruelle, M. Sasaki and Y. Sendouda. larXiv:0803.2742l [gr-qc]. 

[18] S. Mukohyama and L. Randall, Phys. Rev. Lett. 92, 211302 (2004); T. Inagaki, S. Nojiri and 
S. D. Odintsov, JCAP 0506, 010 (2005) j |arXiv:gr-qc/0504054| ; A. D. Dolgov and M. Kawasaki, 
|arXiv:astro-ph/0307442 

[19] O. Bertolami, C. G. Boehmer, T. Harko and F. S. N. Lobo, Phys. Rev. D 75, 104016 (2007); 

T. Koivisto, Class. Quant. Grav. 23, 4289 (2006). 
[20] V. Faraoni. larXiv:0710.T29"T1 [gr-qc]. 

[21] O. Bertolami and J. Paramos. larXiv:0709.3988l [astro-ph], 

[22] S. Nojiri, S. D. Odintsov and P. V. Tretvakov. larXiv:0710.5232l [hep-th]. 

[23] K. Bamba and S. D. Odintsov, JCAP 0804, 024 (2008) |arXiv:0801.0954l [astro-ph]]. 

[24] K. Bamba and M. Sasaki, JCAP 0702, 030 (2007); K. Bamba, ibid. 0710, 015 (2007). 

[25] I. T. Drummond and S. J. Hathrell, Phys. Rev. D 22, 343 (1980). 

[26] A. Raya, J. E. M. Aguilar and M. Bellini, Phys. Lett. B 638, 314 (2006); F. A. Membiela and 

M. Bellini. larXiv:0712.3032l [hep-th], 
[27] L. Campanelli, P. Cea, G. L. Fogli and L. Tedesco. larXiv:0802.2630l [astro-ph], 
[28] M. S. Turner and L. M. Widrow, Phys. Rev. D 37, 2743 (1988). 

[29] K. -T. Kim, P. P. Kronberg, P. E. Dewdney and T. L. Landecker, Astrophys. J. 355, 29 

(1990); K. -T. Kim, P. C. Tribble and P. P. Kronberg, ibid. 379, 80 (1991); T. E. Clarke, 

P. P. Kronberg and H. Boehringer, ibid. 547, Llll (2001). 
[30] Y. Sofue, M. Fujimoto, and R. Wielebinski, Annu. Rev. Astron. Astrophys. 24, 459 (1986); 

P. P. Kronberg, Rept. Prog. Phys. 57, 325 (1994); D. Grasso and H. R. Rubinstein, Phys. Rept. 

348, 163 (2001); A. D. Dolgov, |aTXiv:hep-ph/0110293[ C. L. Carilli and G. B. Taylor, Annu. 

Rev. Astron. Astrophys. 40, 319 (2002); L. M. Widrow, Rev. Mod. Phys. 74, 775 (2003); 

M. Giovannini, Int. J. Mod. Phys. D 13, 391 (2004); V. B. Semikoz and D. D. Sokoloff, ibid. 

14, 1839 (2005); M. Giovannini, |arXiv:astro-ph/0612378~ 
[31] L. H. Ford, Phys. Rev. D 40, 967 (1989). 

[32] A. Golovnev, V. Mukhanov and V. Vanchurin. larXiv:0802.2068l [astro-ph], 

[33] W. Zhao and Y. Zhang, Class. Quant. Grav. 23, 3405 (2006); Phys. Lett. B 640, 69 (2006); 

Y. Zhang, T. Y. Xia and W. Zhao, Class. Quant. Grav. 24, 3309 (2007). 
[34] T. Y. Xia and Y. Zhang, Phys. Lett. B 656, 19 (2007). 

[35] E. Elizalde, J. E. Lidsey, S. Nojiri and S. D. Odintsov, Phys. Lett. B 574, 1 (2003) 

[arXiv:hep-th/0307177] . 
[36] C. Armendariz-Picon, JCAP 0407, 007 (2004). 
[37] H. Wei and R. G. Cai, Phys. Rev. D 73, 083002 (2006). 
[38] H. Wei and R. G. Cai, JCAP 0709, 015 (2007). 
[39] J. B. Jimenez and A. L. Maroto. larXiv:0801. 14861 [astro-ph]. 

[40] T. Koivisto and D. F. Mota, larXiv:0707.0279l [astro-ph]: larXiv: 080 1.36761 [astro-ph]. 
[41] V. I. Demchik and V. V. Skalozub, Eur. Phys. J. C 25, 291 (2002); 

V. V. Skalozub and A. V. Strelchenko, ibid. 33, 105 (2004). 
[42] S. G. Matinyan and G. K. Savvidy, Nucl. Phys. B 134, 539 (1978); H. Pagels and E. Tomboulis, 



20 



ibid. 143, 485 (1978). 
[43] S. L. Adler and T. Piran, Rev. Mod. Phys. 56, 1 (1984). 

[44] M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory (Westviews 

Press, Boulder, Colorado, 1995). 
[45] A. B. Balakin and J. P. S. Lemos, Class. Quant. Grav. 22, 1867 (2005). 
[46] A. B. Balakin, H. Dehnen and A. E. Zayats, Phys. Rev. D 76, 124011 (2007). 
[47] P. Tzanavaris, M. T. Murphy, J. K. Webb, V. V. Flambaum and S. J. Curran, Mon. Not. Roy. 

Astron. Soc. 374, 634 (2007). 
[48] R. A. Battye, R. Crittenden and J. Weller, Phys. Rev. D 63, 043505 (2001). 
[49] P. Stefanescu, New Astron. 12, 635 (2007). 

[50] R. Khatri and B. D. Wandelt, Phys. Rev. Lett. 98, 111301 (2007). 

[51] L. Bergstrom, S. Iguri and H. Rubinstein, Phys. Rev. D 60, 045005 (1999). 

[52] P. P. Avelino et al, Phys. Rev. D 64, 103505 (2001). 

[53] Y. Fujii and M. Sasaki, Phys. Rev. D 75, 064028 (2007). 

[54] E. Garcia-Berro, J. Isern and Y. A. Kubyshin, Astron. Astrophys. Rev. 14, 113 (2007). 
[55] E. W. Kolb and M. S. Turner, The Early Universe (Addison- Wesley, Redwood City, California, 
1990). 

[56] W. L. Freedman et al., Astrophys. J. 553, 47 (2001). 

[57] S. Nojiri and S. D. Odintsov, J. Phys. Conf. Ser. 66, 012005 (2007) |arXiv:hep-th/0611071| . 

[58] K. Dimopoulos and M. Karciauskas. larXiv:0803.3041l [hep-th], 

[59] K. Dimopoulos, Phys. Rev. D 74, 083502 (2006); 76, 063506 (2007). 



21 



